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I. INTRODUCTION 

The calculation of ground state properties of a quan- 
tum system with many degrees of freedom has been ex- 
plored by several means. Exact diagonalisation is usually 
limited to fairly small sized systems. Monte Carlo meth- 
ods are hard to extend to zero temperatures and/or are 
seriously hampered by sign problems in the wave func- 
tion (in particular for fermionic degrees of freedom) . Re- 
cently White has introduced a new algorithm, which 
bears some analogy with the renormalization technique 
in the sense that wave functions of larger systems are 
constructed hierarchically from smaller components. It 
has received the name Density Matrix Renormalization 
Group (DMRG) although the group character is nowhere 
present and even the link with renormalization as induced 
by spatial rescaling is rather weak. The DMRG-method 
has achieved remarkable accuracy for a number of sys- 
tems notably those of a 1-dimensional {d = 1) character. 
Although there is a physically acceptable rational for the 
algorithm, its limitations are not well understood. In par- 
ticular the restriction of the success to d = 1 systems is of 
an empirical nature while theoretically the renormaliza- 
tion idea would equally well work in higher dimensions. 
As noted earlier the method performs poorest |^,^ near 
a quantum phase transition. According to Ostlund and 
Rommer Q this has to do with the hierarchical nature 
of the ground state wave function which is at odds with 
the algebraic correlations in a critical system. 

In this paper we study the Ising model in a Transverse 
Field (ITF) as a model system for a quantum phase tran- 
sition on a 2 dimensional square lattice. It has the advan- 
tage to be exactly soluble in d = 1 dimensions which helps 
checking accuracies for the method we use. Straightfor- 
ward application of the DMRG-methods yields highly ac- 
curate results in this case. The real challenge is d = 2 
dimensions where the phase transition has the same com- 
plexity as that of a, d — 3 dimensional classical Ising 
model. Here straightforward, brute force application of 
the DMRG-technique does not yield convincing results 



and sophistication is called for. Still some results on 
d = 2 dimensional spin system have been achieved by 
White @. 

The new ingredient in this paper is that we grow the 
system with a whole band per step in stead of just one 
site (see figure |l|) . This approach allows us to implement 
the translational symmetry in the VF-direction in a simi- 
lar fashion as Xiang has done for the Hubbard model. 
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FIG. 1. The systems we consider are of dimensions L x W 
where L = 2W, 3W, 5W, 20 and = 2, . . . , 6. The system 
contains three parts: A left-hand part A (shaded), an inter- 
mediate band B (black) and a right-hand part C (white). 
Actually the figure obscures that for L — 2W, JiW, 20 the sys- 
tems are periodical in both directions. At every DMRG step 
part A and B axe contracted. 



The lay-out of this paper is as follows: First we in- 
troduce the ITF. The critical behaviour of this model is 
discussed next. After that, we shift our attention to the 
DMRG. We make a link with perturbation theory and 
show the limitations of the method as a consequence of 
the environment, part C. Afterwards our implementa- 
tion is described. Finally, we present the results of the 
calculations: On one hand, the accuracies achieved and 
on the other hand, the critical properties of the ITF. 
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II. THE ITF 

Since the beginning of the 1960s, the ITF has been 
studied. In first instance, the purpose was to model spe- 
cific materials like KH2PO4 crystals. With the introduc- 
tion of the renornialization group in the 1970s, an other 
interest in this model has arisen. As there exist a simple 
relation to a classical system, the ITF is used as a vehi- 
cle to extend the knowledge of critical phenomena from 
classical spin systems to quantum systems. Chakrabarti, 
Dutta and Sen |^ have recently summarised the prop- 
erties of the ITF. For further details we refer to them. 
We will only mention those properties here, that are of 
explicit use to our calculations. 



A. The Model 

Consider a two dimensional square lattice with length 
L and width W. The lattice is periodic in both directions 
and each lattice site contains a spin-^. The Hamiltonian 
is given by 

L w 

'^ITF = E E i-J^^.,(^t^i,, + + HSl^) (2.1) 

i=i 3=1 

where the 5" are the usual Pauli spin matrices satisfying 



a,/3,7 = x,y,z (2.2) 



As only the ratio H/ J is important, we fix J at J = 1 
and take H > { H < being equivalent). 

This model is translation and reflection symmetric 
in both directions. Moreover, the symmetry operation 
Sfj -Sfj, Sfj ^ -Sf j and S^j leaves the 

model invariant. The operator associated with this is 
S = e-xp{iTT/2{J2i j S^j + LW)). It samples the total 
number of spins pointing upwards and returns whether 
it is odd S = —1, or even S = +1. We call this S the 
spin-reversal operator. The ground state is an eigenfunc- 
tion of all these symmetries operations. 

If — > 0, we end up with a simple 2D Ising model. 
The ground state is degenerate; all spins point either up 
or down in the iS^-direction. The associated phase is the 
classical ordered phase. By taking a rotation in the low- 
est energy space, we can obtain states that are even in 
spin-reversal {S = +1) or odd {S = —1). In the other 
extreme, 1/H —>■ 0, free spins in an external field remain. 
The ground state is unique and has all spins pointing 
down in the iS^-direction. This is the reference state for 
the quantum disordered phase and has value 5 = +1. 
The lowest excitation differs from the ground state by the 
reversal of one spin. So it belongs to the class S = —1. 
We will extensively study the energy gap A between the 
lowest excitation (in S = —1) and the ground state (in 

S = +1); A = i?ex — Egr- 



There is a phase transition between the classical or- 
dered and the quantum disordered state. A clear signa- 
ture of this phase transition is the disappearance of the 
gap A, which occurs for a critical value H = H^- 



B. Critical Behaviour 

As mentioned before, the ITF is closely related to 
a classical Ising model. It can be mapped onto an 
anisotropic Ising model in one dimension higher. In the 
current situation the resulting classical model is of size 
LxW X 00. It contains a weak coupling K^^ in the L x W- 
plane and a strong Ising coupling K]^]^ in the remaining 
direction (exp(— i^Tij) = eH , K± = eJ with e <C 1). 
Chakrabarti et al. [|6| give an overview of the procedure. 
For our purposes the most important consequences are: 

• The correlation length ^ in the strong coupling di- 
rection corresponds to the inverse of our gap A 

(e^A-i). 

• The reduced temperature < = (T — Tc)/Tc corre- 
sponds to our reduced field h = {H — He) /He- 
{t ^ h). 

The well-known relation ^ ^ t~^'^ transforms into A ~ 
h^-'^ . In the 3D anisotropic Ising model the dynamical 
exponent z = 1. We use this in the further discussion. 
The standard finite size scaling methods as described 
in can be applied here. The classical scaling rela- 
tion ^(i, H/-1) = b^ib^/^t^hW-^), for fixed aspect ratio, 
L/W = constant, becomes 



A(/i,VK-i) = \/\{b^'''h,bW- 



(2.3) 



We may set b = W and obtain the scaling expression 

W/^[h, W-^) ^ J^iW^'^h, 1), (2.4) 

showing that W IS. only depends on the combination 
W^^^h. So for /i = all lines W /S. cross at the same 
value; 



1FA(0,I^-^) = A(0,1) 



(2.5) 



This gives us He- If we differentiate ( p.4[ ) with respect to 
h and set h — Q afterwards, we obtain 
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1 - - lOo , 



W \ / A^(0,VK-i) 

= - log ' 



A^(0,(W^ + l)-i) 



(2.6) 



From this we can extract v. 
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III. THE DMRG METHOD 

The DMRG method was first formulated by White 0. 
As it is not a renormahzation group method in the tra- 
ditional sense, it could perhaps better be named an it- 
erative basis truncation method. Gehring, Bursill and 
Xiang 1^ provide an excellent introduction in the appli- 
cation of the DMRG to ID spin systems. Here we will not 
be so extensive. Two important features of the method 
are discussed and our approach is outlined. 



be reproduced. A reformulation of this is: Consider a 
I a) such that {Pj\a) = for all 1(3 j). We know that 
|a) =E,^fK>, thus 



and 



(/3,|a)=^, 



E 



i'j 







v-^ = 0. 



(3.5) 



(3.6) 



A. Limitations by the Environment 

The essence of the DMRG can be described as follows: 
Consider a system consisting of 2 parts, A and B. More- 
over, suppose we have an approximate ground state wave 
function of the combined system |(/)o) = Xlij "^ijWIj)- 
The bases and {|j)} do not have to be complete. 

We want to reduce the number of basis states in part A, 
preserving the ground state wave function |^o) as well as 
possible. |0o) can be expanded as 



'a J I 



where {|a)} spans only a subspace of 
the ground state means that 



(3.1) 



Preserving 



is minimal. 



(3.2) 



The solution to this problem can be obtained by means 
of simple algebra [Q ; Construct the density matrix pui — 
Tlij and select the eigenvectors v°' with the largest 



eigenvalues A" ; p-v" = 
by: \a)=Y.^vf\i). 
A truncation error 



P 



The new basis is now given 



(3.3) 



is introduced to give an indication of the effectiveness of 
the procedure. On basis of experience this truncation er- 
ror is said to be a measure of the error in the calculated 
energy with respect to the exact result. 

There is a peculiarity that was only briefly mentioned 
by White 0]. Suppose we want to use this selection 
scheme to obtain as many states in part A as we already 
have in part B, presuming that there were more states 
in A initially. Define \j3j) = X^i'/'yN)- The ground state 
can be transformed into this set; 



(3.4) 



Thus by orthonormalising the set {|/3j)} we obtain a 
basis set for A in which the wave function can exactly 



j7" is a zero eigenvector of p. Keeping the subspace 
spanned by the \(3j) would make the truncation error p 
equal to zero. This lack of choice only becomes worse 
in case symmetries are implemented; not just the total 
number of non-zero eigenstate is fixed, but even within 
a specific symmetry class the number of non-zero eigen- 
state is dictated by the states in the environment. Later 
on, we will make this explicit for the systems we consider. 



B. The Connection with Perturbation Theory 

Liang and Pang |^ mention that for a given accuracy, 
the number of states needed in a single part of the sys- 
tem grows exponentially with the width of the system. 
At the phase transition, we confirm this observation (fig- 
ure ||). Moreover, both in the weak and strong field limits 
{H <^ 1 and 1/H <C 1), we find a fast convergence, which 
can be explained by perturbation theory. 

Consider the quantum disordered phase. Split the 
Hamiltonian into a unperturbed part Tio = H ^ S^j 
and a perturbation V = - + We 

split the periodical, rectangular system of size L x W 
again in two parts; A and B of sizes I x W and {L — l)xW 
where I is an arbitrary length smaller than L. They both 
contain 2W spins that border the other part. The unper- 
turbed ground state |0) has all spins pointing down in the 
5^-direction. It is the direct product of two equivalent 
states restricted to A and B; |0) = |0)yi|0)B. We know 
that WolO) = -HLW\0) = Ea\0). Perturbation theory 
yields 



\M - |o) + 



(3.7) 



The perturbation flips a pair of neighbouring spins. This 
pair can be in a single part or it can cross the border 
between both parts. In the latter case the spins are adja- 
cent across the boundary between part A and B. Define 
{ I a) a} to be the set of states with the flipped pair in part 
A. Idem for Moreover let {|n)^} be the set with 

one spin flipped on the nth boundary site with B and de- 
fine in an equivalent manner {|n)_B}. The perturbation 
expansion can now be rewritten 
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\ a b n / ^ 

a / \ b / n 



\n)A\n)B + 



(3.8) 



As » 1, it is necessary to reproduce all these terms 
for an accuracy which is equivalent to the first order per- 
turbation theory. The minimal number of states needed 
in part A is therefore 1 for the first term in ( |3.§| ) plus 
2W for all the boundary terms. We have confirmed this 
prediction explicitly in both the small and large H limit 
( H = 1/50,50). 

The same line of reasoning also holds for the second 
and higher order perturbation terms. We expect for an 
error comparable to the nth order perturbation theory 
that m ^ W^", 5E ^ (l/i?)". This is always an upper 
bound for number of states m needed, m < W" for a 
given accuracy 6E ~ (1/7J)". Only when the different 
orders in perturbation theory become distinguishable in 
size - the limit of large H- the equivalence holds. 



C. Exploiting the Symmetries 

We consider systems of sizes L x W. The length L is 
either a muhiple oiW, L ^ 2W, 3VF, 5W, or it is fixed, 
L = 20. The width W is varied from W ^ 2 to W ^ 6. 
The maximal system we study thus contains 6 * 30 = 180 
spins. For L — 2W, 3W and 20 a torus is constructed by 
imposing periodical boundary conditions in both direc- 
tions. For L = 5W the system follows the figure |l| more 
genuinely; it is periodical in the width-direction and open 
in the length-direction. The system is split in a left-hand 
and a right-hand part, both containing m states. A in- 
termediate band, containing the complete basis of 2^ 
states, separates them. This is depicted in figure 0. 

The Hamiltonian of such a system contains many 
symmetries that we can incorporate in our calculation. 
The general form of the included symmetry operators 
is that they are the direct product of three compo- 
nents. Each component acts on one part of the sys- 
tem only. For example, consider the translation oper- 
ator T in the width-direction. This operator is the di- 
rect product of three translations in the individual parts; 
T = TaTbTc- The same holds for the refiection TZ in 
the same direction,?^ — TZaTZbT^c j ^^nd the spin-reversal 
operator S = exp(i7r/2(X]i j >5f ^ -I- LW)) = SaSbSc- 

The ground state |0o) of ths system is translational, 
refiection and spin-reversal invariant; T|0o) = ^|0o) — 
S\(j)o) — \4>o). For systems of infinite size in the clas- 
sical ordered region {L,W oo and H <C 1), it will 
become degenerate with a state that is spin-reversal anti- 



r 



symmetric. In order to take advantage of the symmetries, 
the bases of part A, B and C are chosen to be eigenvec- 
tors of the symmetry operators T and S. TZ is used later 
on. So if {|a)}, {\h)}, {\c)} are the bases of the individ- 
ual p arts ( not to be mixed-up with the notation used in 
(HI)) then 



T^la) = e'''°|a) , SaW) ^ Sa\a). 



(3.9) 



Similar relations hold for the other two sets. Thus 

\<l>o) ^Y.'t'-bc\a)\b)\c) (3.10) 



'ibc 



and appl ication of the symmetry operations together 
with (|J) yields: 



+ fcfc + fee = mod 27r, 

^a'^b'^c — 1- 



(3.11) 
(3.12) 



It is also possible to set up the program to find the lowest 
state in other symmetry classes by forcing other values 
than and 1 in the equations above. 

The Hamiltonian can be written as the sum of Hamil- 
tonians within the separate parts: 'Ha,'Hb and Tic com- 
bined with interactions between parts: TLab,'Hbc and 
Ti-CA', Ti- = Ti-A + T^B + "He + "Hab + Ti-BC + Ti-cA- To 
show how to implement the symmetries, we will discuss 
one element of both types. 

First Ha- It is translational and spin- reversal invari- 
ant, thus 



I'lHAla) = {a'\T^^HArA\a) = 



a'\nA\a) 



{a'lS^^HASAla) = Sa'Sa{a'\TiA\a) 
{a'\HA\a)Ss,,sJk,,ka- 



(3.13) 



It only contains elements within symmetry classes, as one 
would expect. 

Second Hab- Once again, it is translational and spin- 
reversal invariant. Moreover it can be written as 



w 



T^AB — — E '^Cn'^i+l.n 



71=1 
W 

= -Y,{TArBr'+'Sl,Sf+,^,{TATBT-' (3.14) 

n=l 

where I is the length of part A. Sfj_ fiips a spin, so 



SfjS + SSf, = 0. Inserting this and (3.£) in (3.14) gives 
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{a'\{b'\nAB\a)\h) = ~W{a'\Sl,\a){b'\Sf^,^,\h) 



(3.15) 



This substantially reduces the computational effort. Fi- 
nally: the reflection operator TZ is used to make matrix 
elements like {a'\S[-j^\a) real. Naturally we could have 
used this last symmetry TZ more, but it only reduces the 
effort by a factor of 4 while making the program far more 
complex. 



D. The Implementation 

Now we focus on the procedure itself. It is tempting 
to use the ID DMRG method directly: a site is replaced 
by a band. The ground state |^o) of the entire system 
ABC is calculated and the optimal basis for block AB is 
selected through the density matrix. However, one runs 
in severe difhculties as a consequence of the first remark 
on the DMRG. It is instructive to reveal the reason: Us- 
ing the notation above, we define |/3c) = J2ab'f'a.bc\o)\b)- 
We know that T|0o) = S\4>q) = {(po), thus 



SaSbIPc 



(3.16) 



The distribution over the symmetry classes in part 
C forces the selected states in block AB to be in 
"conjugate" -classes. To overcome this problem, we need 
to increase the number of states in part C. In that case 
we can really make a selection and shift into important 
symmetry classes. 

In the ID procedure the solution is to add one extra 
site to the environment. The number of states in the 
environment is then doubled. In our set-up this would 
correspond to adding an extra band between B and C. 
This is computational far too expensive. We now in- 
troduce variants on White's infinite-size and finite-size 
algorithms that increase the number of states in the 
part C. 

First we consider our infinite size approach. We only 
have to describe one step in the process as it is an induc- 
tive method. We have a basis of m states for a system of 
length I. 

• We construct the combined system as depicted in 
figure ||-a by taking this basis in part A and C to- 
gether with the complete basis in the intermediate 
band B. {L = 21 + 1) 

• We calculate the ground state \<j>o) and obtain m 
basis states for a system of length ^ -I- 1 by orthonor- 
malising {\(3c)}. 

• Suppose that block AB has / symmetry classes. To 
every symmetry class we add m/ f basis states con- 
structed randomly from the m2'^ states in A and 



B. We end up with m + f ■ m/ f 
for a system of length ^ -I- 1 . 



2m basis states 



• In part A wc now take the m basis states for a sys- 
tem of length I and in part C we take the newly 
constructed 2m states for length l + l. (L ~ 21 + 2) 
This yields the configuration in figure ^b. 

• We calculate the ground state \(j)o) and obtain 2m 
basis states for length I + 1 hy orthonormalising 
{|/3c)}. We replace the basis of part C by this basis 
and repeat this step a couple of times (^ 3). 



• We select from the 2m basis states for length I 
m states on basis of the density matrix. 



1 



Now we have returned to the original situation with the 
exception that I has increased by one. The new ingre- 
dient is thus to add m random states to the basis and 
iterate until the result has converged. 



a) start-up 



b) iteration 
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I 1 l+l 

FIG. 2. A inductive step in the infinite-size procedure con- 
sists of a start-up to obtain an initial approximation for states 
in a system of length I + 1 and iterative calculations to make 
the basis converge. The numbers in the rectangle are the num- 
ber of states in the parts. The intermediate band B always 
contains the complete basis of 2^ states. 



In the same line our finite size approach lies. Suppose 
we have basis sets of m states for lengths l,L — I — 1 and 
L — Z — 2, where L is now fixed and independent of I. 

• We take the basis for I in part A, the basis for 
L — I — 1 in part C and the complete basis of the 
band in part B. See figure ||-a. 

• We calculate the ground state |0o) and obtain a 
basis for length Z -I- 1 by orthonormalising {|/3c)}. 

• In the same way as in the infinite-size algorithm we 
add m randomly chosen states to this basis. 

• In part C we take the 2m basis states for length 
I + 1 and in part A the m states for length L — l — 2. 
This is depicted in the first of the two pictures in 
figure p|-b. 
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• We calculate the ground state \4>o) and obtain 2m 
basis states for L — Z — 1. 

• In part C we take the 2m basis states for length 
L — I — 1 and in part A the m states for length 
see the second picture in figure |^b. 

• We calculate the ground state \(j)o) and obtain 2m 
basis states for I + 1. These last four steps are re- 
peated a couple of times (~ 3) 

• We select from the 2m basis states for length / + 1 
m states on basis of the density matrix. 

Once again we have returned to our starting position 
while increasing the length I by one. By sweeping 
through the system we can therefore systematically im- 
prove the basis. This method convergences at a similar 
speed as the ID approach; After 3 sweeps through the 
system the final result is achieved. 
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a) start-up 
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2m 





I 1 L-l-1 

FIG. 3. A inductive step in the finite-size procedure also 
consists of a start-up to obtain an initial approximation for 
states in a system of length / + 1. Afterwards we move back 
and forth between lengths I and Z + 1 to make this converge. 

The computational effort scales as m^L^. In gen- 
eral m > 2^. This clarifies the bound on the width. 
The alternative is to follow Liang, Pang |^ and White 
by adding one site per step. We can then still use S. 
The calculation would scale as m^L^r- Our approach 
includes the symmetry requirements of the ground state 
and up to 8 it is similar in speed as theirs. Applying 
our method to models where the number of particles or 
total spin is conserved instead of S, the calculations can 
been substantially reduced and systems of widths W > 6 
are pulled within reach. 

The largest calculation, W = 6, L = 18, m = 200 took 
48h of computer time per 7?- value (at 462 SPECfp92). 
To determine the gap A two such points are needed 
{S = +1 and S = -1). 



IV. RESULTS 



We have performed two kinds of calculations: First, we 
made a check on the accuracy of the method. Second, we 
have calculated the gap A for various widths W, aspect 
ratio's L = xW and fields H in order to find through 
finite size scaling the phase transition point He and the 
critical exponent ly. 



A. The Accuracy of the Energies. 



The strict method to determine the error in the energy 
5Em for given number of states m is to compare the re- 
sults with the exact value Eg^; SE„i = E,n — Eg^. 
This would limit us to small systems of sizes comparable 
to 6 X 6. In the literature |^ it is noted that the error 
6E decreases exponentially with the number of states m 
included. We confirm that statement explicitly for these 
small systems. Moreover we use this feature to test the 
accuracy for far larger systems. The energy Em is com- 
pared with the result for a larger number of states. For 
instance m < 128; SEm ~ Em — i?i28- The error SEm is 
largest near the phase transition as can be seen in figure 




a 



o 

; -5.0 - 



-6.0 - 



-7.0 ^ ' ' ' ' ' ' ' 1 

1.0 2.0 3.0 4.0 5.0 

H 

FIG. 4. The accuracy of the DMRG method for different 
of states m (numbers in graph) as function of the field H. 
The system is periodical in both directions with dimensions 
W — 4: and L — 20. The reference value is taken from a 
DMRG calculation with m = 64. 
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As the phase transition occurs near iJ = 3, we take 
H — 3 a,s a.n example to study the dependence of the 
error SE on the width W. The error dEm increases ex- 
ponentially with growing width W (figure ^ . 

-2.0 I ; ; ; ; ; 1 
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-7.0 
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x16 
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x64 



x16 



x32 



X64 



x32 



X64 



x64 



x32 



4 
W 



FIG. 5. The accuracy of the DMRG method for given num- 
ber of states m (numbers in graph) as function of the width 
W. H — 3 and L = 20. The system is periodical in both 
directions. The reference value is taken from a DMRG calcu- 
lation with m — 128. 




3.10 



FIG. 6. The scaled gap WA as function of the field H for 
aspect ratio L = 2W. W = 2,3, 4, 5, 6. The curves become 
steeper with increasing width W. The crossings for consec- 
utive widths are encircled. The system is periodical in both 
directions, m = 128, for W = 2,3,4,5 and m = 200 for 
W = 6) 



B. The Phase Transition and the Critical Exponent. 



The ph ase transition point He is determined through 
equation ( ^.5[ ). We plot WA versus H (figure || and Q)- 
The curves would intersect precisely at He , if it were not 
for corrections to scaling. These become quit e large when 
W = 2,3. Afterwards we use formula (2^) to obtain ly 
at the intersection of the curves for consecutive widths 
W. The results are listed in table |. For W = 6 and 
L — 2W, 3W we are at the limit of our precision, when 
we take m — 128 states. We therefore set m — 200 in 
this case. 

Apart from these periodical systems, we have also con- 
sidered systems where the periodical connection between 
A and C is removed. The removal of this boundary con- 
nection has two effects: First, the accuracy of the cal- 
culated energies will increase substantially as the size of 
the interacting boundary is halved. Second, the correc- 
tions to scaling will increase. To make up for this second 
effect, we have to resort to fairly large systems; L = 5W. 
This is depicted in figure | . 

From the values in table | we note that the corrections 
to scaling for v are still fairly large for these system sizes 
5%). We found that these corrections could not be 
compensated by introducing an irrelevant scaling field in 
the relation (O). 
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FIG. 7. Idem as figure H with now the ratio L = 3W. 
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FIG. 8. The scaled gap AW for systems with open bound- 
ary conditions in the length-drections and periodical in the 
width-direction. L = 5W and m = 64. 





L = 


2W 


L = 


■3W 


L = 


5W 


w 


He 


V 


He 


V 


He 


V 


2-3 


3.113 


0.74 


3.110 


0.73 


3.101 


0.74 


3-4 


3.068 


0.69 


3.067 


0.68 


3.062 


0.69 


4-5 


3.054 


0.67 


3.053 


0.67 


3.051 


0.67 


5-6 


3.049 


0.66 


3.047 


0.65 


3.046 


0.66 



TABLE I. 



and the crit- 



The phase transition point He 
ical exponent v. We take He to be the value where 
WA{W'^) = {W + 1)A{{W +1)-^). V is calculated through 
equation (2.6). The first two aspect ratio's L — 2W,3W are 
with periodical boundary conditions, the last L — 5W is with 
open boundary conditions in the length-direction. 



V. CONCLUSION 

In this paper we have presented an adaption of the 
DMRG method to two dimensional spins systems. We 
follow the route of adding complete bands instead of sin- 
gle sites to the system. The latter was done by Liang, 
Pang and White This modification allows us to 
force a translational symmetry in the width-direction. 
The advantage of implementing this symmetry is that a 
ground state with specific translational properties can be 
targeted. Moreover, the space in which the ground state 
is sought is reduced substantially. This is especially use- 
ful in systems with Goldstone modes or similar gapless 
excitation spectra where the lowest excitations belong to 
different symmetry classes than the ground state. 

The computational effort still remains similar to the 
approach of adding single sites as the larger space of the 
band (2^ instead of 2) is offset by three reductions: First, 
the ground state can be written more compactly (a factor 
W reduction). Second, we only need to apply one opera- 
tor per boundary instead of W operators. Third, the 
sub-system (part A) grows with a full band instead of a 
single site per step (factor W). 

We have only considered systems of widths upto W = 
6. In models where the total spin or the number of par- 
ticles is conserved, we can go to larger widths. 

We observe that at criticality, the number of states m 
needed for a given accuracy SE/E grows exponentially 
with the width W, in full agreement with Liang and Pang 
. Moreover, we have proven that far enough from the 
phase transition the method will reproduce perturbation 
theory. 

The procedure does not get stuck at local minima as 
was sometimes experienced by White and Scalapino [ p^ . 

The gap A we have calculated is a nice example of 
the use of symmetry classes. The results for the critical 
properties, He = 3.046 and v = 0.66, are in reasonable 
agreement with the series expansions of Pfeuty and El- 
liott and with the cluster Monte Carlo calculations 



of Blote [H. 

As yet, this method is not as accurate as the more 
traditional methods like Monte Carlo simulations. The 
accuracy could be improved when a larger width could 
be handled by including several hunderds of states. At 
present this would require the use of a super-computer. 
Still it has to be stressed that the DMRG can handle 
problems that are out of reach of Monte Carlo simula- 
tions due to the "sign-problem" . 
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